Crustal structure and apparent tectonic underplating from receiver function analysis in South Island, New Zealand

We utilize seismic converted phases on more than 700 receiver functions calculated for 42 stations in the South Island, New Zealand, to infer crustal and uppermost mantle structure. We determine the crustal thickness from direct observations of conversion from the Moho interface and infer zone of the maximum thickness being located along the axis of the Southern Alps, just east from the Alpine fault. The crustal root widens from north to south in the direction perpendicular to the Alpine fault and appears to have an asymmetric structure. Stations in the alpine portion of island show evidence for prominent midcrustal conversions. Significant crustal thickening is developed in response to both the convergent component of the motion on the Alpine fault and subduction in the Fiordland region. We propose two models for a strong uppermost mantle conversion that occurs at depths between 33 and 83 km on 16 stations and forms a large continuous feature along the east coast and in the central portions of the South Island. Our preferred model attributes upper mantle conversion to tectonically underplated oceanic crust formed by late Oligocene-Miocene spreading between the Australian and Pacific plates, which was detached from the Australian plate and tectonically underplated under the South Island. An alternative model attributes the upper mantle conversions to long-lived seismic fabric created by subduction of the Gondwanaland margin

P wave velocity variations in the Coso Region, California, derived from local earthquake travel times

Inversion of 4036 P wave travel time residuals from 429 local earthquakes using a tomographic scheme provides information about three-dimensional upper crustal velocity variations in the Indian Wells Valley-Coso region of southeastern California. The residuals are calculated relative to a Coso-specific velocity model, corrected for station elevation, weighted, and back-projected along their ray paths through models defined with layers of blocks. Slowness variations in the surface layer reflect local geology, including slow velocities for the sedimentary basins of Indian Wells and Rose valleys and relatively fast velocities for the Sierra Nevada and Argus Mountains. In the depth range of 3–5 km the inversion images an area of reduced compressional velocity in western and northern Indian Wells Valley but finds no major velocity variations beneath the Coso volcanic field to the north. These results are consistent with a recent study of anomalous shear wave attenuation in the Coso region. Between 5 and 10 km depth, low-velocity areas (7% slow) appear at the southern end of the Coso volcanics, reaching east to the Coso Basin. Numerical tests of the inversion's resolution and sensitivity to noise indicate that these major anomalies are significant and well-resolved, while other apparent velocity variations in poorly sampled areas are probably artifacts. The seismic data alone are not sufficient to uniquely characterize the physical state of these low-velocity regions. Because of the Coso region's history of Pleistocene bimodal volcanism, high heat flow, geothermal activity, geodetic deformation, and seismic activity, one possibility is to link the zones of decreased P velocity to contemporary magmatic activity

A Restoration Method for Impulsive Functions

A method is presented for enhancing the resolution of impulsive functions which have been degraded by a known convolutional disturbance and by the addition of white noise. An autoregressive model is employed to represent the spectrum of the ideally resolved impulsive function. The method is flexible in that it allows constraints to be incorporated into the resolution scheme. Two quite diverse examples are presented as illustration

Crustal structure of the Borderland-Continent Transition Zone of southern California adjacent to Los Angeles

We use data from the onshore-offshore component of Los Angeles Region Seismic Experiment (LARSE) to model the broad-scale features of the midcrust to upper mantle beneath a north-south transect that spans the continental borderland in the Los Angeles, California, region. We have developed an analysis method for wide-angle seismic data that consists primarily of refractions, lacks near-offset recordings, and contains wide gaps in coverage. Although the data restrict the analysis to the modeling of broad-scale structure, the technique allows one to explore the limits of the data and determine the resolving power of the data set. The resulting composite velocity model constrains the crustal thickness and location and width of the continent-Borderland transition zone. We find that the mid to lower crust layer velocities of the Inner Borderland are slightly lower than the corresponding layers in the average southern California crust model, while the upper mantle velocity is significantly higher. The data require the Moho to deepen significantly to the north. We constrain the transition zone to initiate between the offshore slope and the southwest Los Angeles Basin. If the Inner Borderland crust is 22 km thick, then the transition zone is constrained to initiate within a 2 km wide region beneath the southwest Los Angeles Basin, and have a width of 20–25 km. The strong, coherent, and continuous Pn phase suggests the Moho is coherent and laterally continuous beneath the Inner Borderland and transition zone. The Inner California Borderland seems to be modified and thickened oceanic crust, with the oceanic upper mantle intact beneath it

An Evaluation of Proposed Mechanisms of Slab Flattening in Central Mexico

Central Mexico is the site of an enigmatic zone of flat subduction. The general geometry of the subducting slab has been known for some time and is characterized by a horizontal zone bounded on either side by two moderately dipping sections. We systematically evaluate proposed hypotheses for shallow subduction in Mexico based on the spatial and temporal evidence, and we find no simple or obvious explanation for the shallow subduction in Mexico. We are unable to locate an oceanic lithosphere impactor, or the conjugate of an impactor, that is most often called upon to explain shallow subduction zones as in South America, Japan, and Laramide deformation in the US. The only bathymetric feature that is of the right age and in the correct position on the conjugate plate is a set of unnamed seamounts that are too small to have a significant effect on the buoyancy of the slab. The only candidate that we cannot dismiss is a change in the dynamics of subduction through a change in wedge viscosity, possibly caused by water brought in by the slab

Modeling path effects in three-dimensional basin structures

Path effects for seismic wave propagation within three-dimensional (3-D) basin structures are analyzed using a reciprocal source experiment. In this experiment, a numerical simulation is performed in which a point source is excited at a given location and then the wave field is propagated and recorded throughout a 3-D grid of points. Using the principle of reciprocity, source and receiver locations are reversed. This allows the modeling of path effects into a particular observation site for all possible source locations using only one simulation. The numerical technique is based on the use of paraxial extrapolators and currently tracks only acoustic waves. However, the method is capable of handling arbitrary media variations; thus, effects due to focusing, diffraction, and the generation of multiple reflections and refractions are modeled quite well. The application of this technique to model path effects for local earthquakes recorded at stations in the Los Angeles area of southern California indicates the strong influence of the 3-D crustal basins of this region on the propagation of seismic energy. The modeling results show that the Los Angeles, San Fernando, and San Gabriel basins create strong patterns of focusing and defocusing for paths into these stations from various source locations. These simulations correlate well with earthquake data recorded at both stations. By comparing these calculations with earthquake data, we can begin to evaluate the importance of these basin effects on observed patterns of strong ground motions

A detailed image of the continent-borderland transition beneath Long Beach, California

New crustal images beneath Long Beach, California show the region of the Inner Borderland to continent transition. The cross-sections are obtained from stacked autocorrelations of virtual sources generated from oil-industry data recorded in the city of Long Beach, CA. They show that the Moho is dipping at 65° and obliquely truncates an ∼10 km thick flat-lying lower crustal fabric. The Moho appears to be fault controlled and an integral part of the extrusion of the Catalina Schist that underlays the Inner Borderland. The basement interface has significant offsets of up to 2 km, none of which correspond to the mapped trace of the Newport–Inglewood Fault

A detailed image of the continent-borderland transition beneath Long Beach, California

New crustal images beneath Long Beach, California show the region of the Inner Borderland to continent transition. The cross-sections are obtained from stacked autocorrelations of virtual sources generated from oil-industry data recorded in the city of Long Beach, CA. They show that the Moho is dipping at 65° and obliquely truncates an ∼10 km thick flat-lying lower crustal fabric. The Moho appears to be fault controlled and an integral part of the extrusion of the Catalina Schist that underlays the Inner Borderland. The basement interface has significant offsets of up to 2 km, none of which correspond to the mapped trace of the Newport–Inglewood Fault

Modeling of residual spheres for subduction zone earthquakes: 1. Apparent slab penetration signatures in the NW Pacific caused by deep diffuse mantle anomalies

We have computed focal residual spheres for 145 subduction zone earthquakes along the northwest edge of the Pacific using regional and global mantle velocity models from tomographic inversions. The mantle models explain much of the observed residual sphere data and, to a certain extent, suggest the location of mantle velocity heterogeneities which are responsible for various residual sphere patterns. For most deep events considered, the fast slablike residual sphere anomalies are caused by diffuse heterogeneities, mainly of deep lower mantle and receiver mantle origin rather than by an extension of the slab. The region immediately below the deepest earthquakes, depths of 650–1500 km, has an effect usually smaller than or comparable to the effect of other regions of the mantle. Without a proper account of the teleseismic effect, attributing the long-wavelength anomalies of the residual sphere to near-source slab effects alone, or even primarily, is not valid. The fast bands in many observed residual spheres agree with seismicity trends. Once the deep mantle and receiver mantle effects are removed, these may give the approximate orientation, but not the depth extent, of near-source fast velocities. For most deep earthquakes under Japan the predominant fast band is subhorizontal rather than near vertical. This type feature would be overlooked in conventional residual sphere studies using only steeply diving rays and cosine weighting of the data

A Born-WKBJ inversion method for acoustic reflection data

Density and bulk modulus variations in an acoustic earth are separately recoverable from standard reflection surveys by utilizing the amplitude-versus-offset information present in the observed wave fields. Both earth structure and a variable background velocity can be accounted for by combining the Born and WKBJ approximations, in a "before stack" migration with two output sections, one for density variations and the other for bulk modulus variations. For the inversion, the medium is considered to be composed of a known low-spatial frequency variation (the background) plus an unknown high-spatial frequency variation in bulk modulus and density (the reflectivity). The division between the background and the reflectivity depends upon the frequency content of the source. For constant background parameters, computations are done in the Fourier domain, where the first part of the algorithm includes a frequency shift identical to that in an F-K migration. The modulus and density variations are then determined by observing in a least-squares sense amplitude versus offset wavenumber. For a spatially variable background, WKBJ Green's operators that model the direct wave in a medium with a smoothly varying background are used. A downward continuation with these operators removes the effects of variable velocity from the problem, and, consequently, the remainder of the inversion essentially proceeds as if the background were constant. If the background is strictly depth dependent, the inversion can be expressed in closed form. The method neglects multiples and surface waves and it is restricted to precritical reflections. Density is distinguishable from bulk modulus only if a sufficient range of precritical incident angles is present in the data